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I. INTRODUCTION 

Catalytically activated processes play a significant role in numerous technologies as they 
serve to produce required products from the species which are non-reactive in normal physical 
conditions; these chemically stable species may, however, enter into reaction in the presence 
of some third substance - the catalytic substrate |I|-£§. Despite of the widespread use of 
such processes, the knowledge of the underlying physics and chemistry still rests largely on 
phenomenological ideas and prescriptions, and thus remains a challenge for fundamental 
research [l]-[5|]. 
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At the simplest, mean-field level of description, reaction kinetics theory presumes that 
the reaction rate should be considered as the product of the reactant densities and the 
rate constant, which is dependent on the nature of the binding forces and on the particles' 
dynamics. This rate constant is proportional to the hopping rate if the process is diffusion 
controlled, or to the reaction frequency, in case when the reaction is kinetically controlled 
f^H- Hence, a clear cut separation is presumed to exist between the local variables, that 
could be derived, say, from quantum mechanics, and collective variables, expressed in the 
most simple way as the product of mean densities of the particles involved. 

Recently, one of the most active developments in the field has been the recognition 
of the substantial importance of the multi-particle effects, spatial fluctuations and self- 
organization, as opposed to conventional local energetic considerations, which emphasized 
the purely chemical constraints and focused on the refined descriptions of the elementary 
reaction act. Subsequently, statistical physics concepts were introduced in order to describe 



anomalous fluctuation- induced behaviors of non-catalytic chemical reactions [ 10 — 25 1 and 



simplest catalytic schemes as well as such to gain a better understanding of such 

collective phenomena as wave formation, presence of excitable media or stochastic aggrega- 
tion in chemical systems (see Refs. [ P8| , |52"U53| and references therein). One conclusion that 
can be drawn from the statistical physics approach is that fluctuations, either spatial or 
temporal, may drive the reactive system into a set of new states, that cannot be understood 
and described in terms of mean-field kinetic equations. 

In this paper we discuss in detail the influence of spatial fluctuations, statistical self- 
organization effects and effects of random diffusive motion of reactants on the kinetics of 
catalytic reactions, using as a particular example the Langmuir-Hinshelwood type reaction 
scheme. This reaction process, which is also often referred to as the monomer- monomer 
catalytic scheme, involves two different kinds of species, A and B, which are deposited 
(continuously in time with mean intensity /) onto the catalytic substrate by an external 
source, diffuse on the substrate and react at encounters forming an inert reaction product 
O, A + B — > O, which is then removed from the system. Our aims here are twofold. First, 



we show that this seemingly simple catalytic reaction (which has, in fact, several practical 
applications (see, e.g. Ref. p6| , p7| )) shows quite a rich behavior and represents an ideal 
illustrative example of the statistical physics effects in the reaction kinetics, which may be 
generic to more complicated schemes involved in real "stuff" catalytic processes. Second, we 
develop a unified analytical description which takes into account explicitly the influence of 
spatial correlations, dimensionality of space and the way how the particles are introduced into 
the reaction bath on the time evolution of observables - mean particle densities. Moreover, 
our approach can be routinely generalized for kinetic description of more complex reaction 
schemes. 

We focus here on several different aspects of the Langmuir-Hinshelwood reaction scheme. 
In particular, we address the question of how the kinetics depends on the dimensionality d 
of the reactive system. In addition to the standard Langmuir-Hinshelwood model in which 
the catalytic substrate is a two dimensional flat surface, we analyse kinetics of A + B — > 
reactions followed by an external input of reactive species in one- and three-dimensional 
systems, i.e. the situations appropriate to reactions in capillary geometries |37[] and annealing 



of the radiation damage in solids [H]. We show that the monomer-monomer catalytic 
reaction proceeds quite differently in low dimensional (d = 1, 2) and three-dimensional 
systems. Besides, we examine how the way of particles' injection into the system affects 
the properties of stationary states and also how these stationary states are approached in 
time. We consider here two different types of external input; in the first one (case I) the A 
and B particles are introduced independently of each other at random moments of time and 
at random positions in space, while in the second case (case II) an A and a B species are 
introduced in correlated pairs of a fixed radius, the pairs being injected at random moments 
of time and at random positions in space. We show that also the way of input does matter 
significantly and may result in a completely different behavior. 

The paper is structured as follows: In Section II we formulate the model, introduce 
a general analytic approach and, in terms of this approach, derive closed-form equations 
describing the time evolution of mean particles' densities and pairwise correlation functions. 



In Section III we present solution to these equations in one-, two- and three-dimensional 
systems corresponding to different ways of particles injection to the reaction bath. Finally, 
we conclude in Section IV with a brief summary of our results and discussion. 



II. DEFINITIONS AND BASIC EQUATIONS. 

In this section we present kinetic description of the monomer-monomer catalytic scheme 
involving diffusive particles in terms of a certain analytical approach, which takes explicitly 
into account the influence of pairwise correlations on the time evolution of mean densities; 
such a description has been first proposed in Ref. ||11|| , which analysed the effects of fluc- 
tuations on the kinetics of A + B — > O reactions and yielded the celebrated t~ d l A -\aw for 
the decrease of particles mean densities. The subsequent works P^j2^j25|j38| -^6| extended 



this approach to more general reaction schemes (e.g. reversible and coagulation reactions), 
three-body and catalytic reactions, and also included the possibility of mutual long-ranged 
interparticle interactions. Here we will focus on the application of this approach to the 
analytical description of the A + B — > O reactions kinetics in systems with a continuous, 
random external input of the reactive species. We will recover some of already known results 
obtained for the monomer-monomer catalytic scheme involving diffusive monomers and will 
establish several new ones concerning the long-time relaxation of mean particles density to 
their steady-state values and the dependence of this steady-state densities on the system 
parameters. 

We start with the formulation of the model. Consider a d- dimensional reaction bath 
of volume V (we suppose that V is sufficiently large such that we can discard different 
finite-size effects, e.g. hard-core exclusion between particles, the termination of reaction 
within a finite time interval or saturation (poisoning)) in which A and B particles are 
continiously introduced by an external random source. The statistical properties of the 
source will be defined below. After injection to the system, the A and B particles begin 
to diffuse. For simplicitly we assume that their diffusion constants are equal to each other, 
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i.e. Da = Db = D. It will be made clear below that such a description is also appropriate 
to the case of non-equal diffusivities; the calculations in this case will be only essentially 
more lengthy [P?| . Now, the reaction event is defined as follows: When any two A and 
B particles approach each other at a fixed separation R (the reaction radius), they may 
enter into reaction forming (an inert with respect to the reaction) reaction product O. 
The recombination upon an encounter of A and B happens with a finite probability p 
(with probability q = 1 — p the particles can be reflected) which defines the constant of an 
elementary reaction act, K. This constant describes the intrinsic chemical activities of A and 
B molecules and is dependent on the nature of the intra- and intermolecular binding forces. 
In the following we will suppose that this purely "chemical" constant is known a priori. 
Besides, we will assume that K is large (i.e. the probability of particles reflection from 
each other is low) and thus will emphasize the "statistical physics" effects on the reaction 
kinetics, rather than the effects of chemical constraints. Consequently, all the factors which 
are exponentially small with K will be neglected here. 

Now, we define the statistical properties of the particles generation more precisely. Let 
J^(r, t) and Js(r, t) be the local, at the point with the radius-vector r, intensities of the 
production rates of A and B particles. We take that the volume average values of the 
production rates obey 

i / drI A (r,t) = 1 / drI B (r,t) = I, (1) 
V Jv V Jv 

which means that particles generation is steady in time and the mean production rates / of 
A and B particles are equal to each other. 

Defining the correlations in the production rates Ia,b{t, t) we will consider here two dif- 
ferent situations. In the first one (case I), we suppose that A and B particles are introduced 
into the reaction bath statistically independent of each other [|16|-|22|| ; that is, the fluctuations 
of the sources are correlated neither in space, nor in time. In this case, we have 

\- I drI A (r,t)I A (r + X,t + r) - I 2 = I 8(X) 5(r), (a) 
V Jv 
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i / drI B (r,t)I B (r + X,t + r) - I 2 = I 8(X) 5(r), (b) 
V Jv 

i J v drI B {r,t)I A {r + X,t + r) = (c) (2) 

In the second case we suppose that A and B particles are introduced as correlated A — B 
pairsQ, separated by a fixed distance X g |T7| , |20| -p2| . Such type of external generation may arise 
in chemical systems in which complex reaction product O is continuously forced to break- 
up by an external radiation (say, laser pulses) into the correlated pairs of the component 
molecules. Here the radius of pair X g will be mainly determined by the difference of energy 
"attributed" to O and the energy required to dissociate the reaction product. Another 
example in which such pairs are produced is the annealing of radiation damage in solids. 
When the solid is irradiated, atoms are knocked out of their places in the lattice to become 
interstitials and leave behind a vacancy, and then the vacancies and interstitials diffuse and 
recombine. In case II, different A — B pairs are statistically uncorrelated and occur, with 
an average intensity J, at random positions in the reaction bath. Then, the fluctuations of 
the sources obey ||T7| , p9| : 



l - [ drI A (r,t)I A (r + X,t + r) - I 2 = I S(X) S(r), (a) 

V Jv 

l - [ drI B (r,t)I B (r + X,t + r) - I 2 = I 5(X) 5(t), (b) 

V jv 

l - J drI B (r,t)I A (r + X,t + r) = ld (X g ) I 5(\X\ - X g ) (c) (3) 



In Eq.(3.c), the parameter 7d(A 9 ) is the normalization factor which arises because of different 
possible angular orientations of a given A — B pair in a (^-dimensional continuum; the value 
of 7d(A s ) depends on the dimensionality of the reaction bath and for d = 1, 2 and 3 equals 
respectively 1, (27rA g ) 1 and (47rA^) _1 . 

Let C A (r, t) and C B (r,t) denote the local densities of A and B particles at point with 
radius- vector r at time t. The time evolution of local densities due to the diffusion of species, 



The particles in the pair can, of course, diffuse apart after injection. 
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their reaction and an external production can be described by the following reaction-diffusion 
equations [p^, p^ , p0| 1 : 

C A {r,t) = - ld {R)K f dr'6(\r-r'\-R)C A (r,t)C B (r',t) + 

+ DA r C A (r,t) + I A (r,t), (4) 



C B (r,t) = - ld {R)K / dr' 5(\r-r'\-R)C B (r,t)C A (r',t) + 

Jv 

+ DA r C B (r,t) + I B (r,t), (5) 

where the symbol A r denotes the (i-dimensional Laplace operator acting on the spatial 
variable r, and the integration with the delta-function S(\r — r'\ —R) accounts for all possible 
orientations of an A — B pair, at which an elementary reaction act can take place. 

Now, an experimentally accessible property is not, however, the local density, but rather 
its volume averaged value 

C(t) = 1 J v drC AiB (r,t) (6) 

To find an equation which governs the time evolution of C(t), let us first represent the local 
densities in the form 

C A , B (r,t) = C(t) + 5C A , B (r,t), (7) 

where SC A)B (r,t) will denote local deviations of particles' densities from their mean values. 
By definition, 

i J v dr5C AB (r,t) = 0. (8) 

Then, substituting Eq.(7) into Eqs.(4) and (5) and taking the volume average, we get the 
following equation 

C(t) = -K [C 2 (t) + G AB {\\\ = R,t)] + J, (9) 
in which G AB (X,t) stands for the pairwise, central correlation function of the form 



Gab{\ t) 




(10) 



the variable A being a d- dimensional correlation parameter. 

Therefore, Eq.(9) shows that the time evolution of the mean particle density is ostensibly 
coupled to the evolution of the pairwise correlations in the reactive system. Neglecting these 
correlations, i.e. setting Gab{\ t) = 0, which is equivalent to the traditional, "mean-field" 
assumption that the particles' spatial distribution is uniform, we obtain the customary, text- 
book " law of mass action" . Such an approximation predicts a linear in time growth of mean 
density at relatively short times, i.e. 



and, in the large-t limit, an exponentially fast relaxation to the equilibrium density C(t — 



The short-time behavior as in Eq.(ll) is quite reasonable and describes the regime in 
which the particles are merely added into the (initially empty) system by the external source 
and the reaction between them is negligible, i.e. the regime in which particle density re- 
mains very small. As for the analytical prediction in Eq.(12), one may question its validity 
on intuitive grounds. First, in the system under consideration diffusive motion of particle 
is the only mechanism to bring particles together and let them react. This motion is es- 
sentially (/-dependent and should evidently entail d- dimensional behaviors, whilst Eq.(12) 
is independent of the dimensionality of the reaction bath. Second, it shows that with an 
increase of the chemical reaction constant K the equilibrium density tends to zero, which 
is apparently an artificial behavior. Below we will show that the actual behavior of C(t) 
as t — > oo is very different. This turns out to depend essentially on the dimensionality of 
the reactive system and also on the way how the particles are injected into the system. We 
show analytically that in the case I, when A and B particles are introduced into the system 



C(t) oc It, 



(11) 



oo) = (I/K) 1 / 2 , i.e. 



C(t) oc (h l/2 \ 1 - exp(-2(/K) 1 / 2 t) + ... 



(12) 
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statistically independ of each other, in low dimensional systems (i.e. d = 1,2) there occurs 
a strong reaction-induced spatial organization of particles. Here the absolute value of the 
pairwise correlation function GAB{R,t) (which is negative) grows in time, which induces 
an unlimited (in absence of hard-core exclusion) growth of mean particles' densities C(t). 
Therefore, in the case I in low dimensional systems the steady state does not exist, in a 
striking contradiction to the prediction of Eq.(12). Similar results have been obtained previ- 



ously in Refs. |16|-P^,p^] . In the case I, in three dimensional systems the steady state density 
is well-defined but, however, turns out to be different from that predicted by Eq.(12). In 
particular, it is dependent on particles' diffusion constant, which reflects the pathological 



behavior of correlations in the reactive system ||T6|J17| , |22|j39| . Actually, we show that correla- 
tions of fluctuations in particles local densities are essentially long-ranged and obey, in the 
steady-state, an algebraic law G AB (\X\,t) ~ 1 / 1 A | . Finally, we proceed to show that in three 
dimensional systems with uncorrelated generation of particles the long-time relaxation of 
mean particles density to its steady-state value is described by a power-law with character- 
istic exponent —1/2, in contrast to the exponential dependence predicted by Eq.(12). Such 
a behavior has been first conjectured in Ref. [^T| and here will be deduced analytically. For 
the case II, we show that steady-state density exists in all dimensions [ |TT| , [2"2"PU| . However, 
its value is also different from that predicted by Eq.(12) and depends, for instance, on the 
radius of the generated pairs, X g . The approach to the steady state in this case is also 
not exponential in time and is described by a universal algebraic law, reminiscent of the 



long-time approach to the equilibrium in reversible chemical reactions 12.13.38-42 



We turn back to Eq.(9) and continue our analysis of the binary reaction kinetics taking 
into account the influence of pairwise correlations on the time evolution of C(t). Omitting 



the details of the derivation, which can be found in Refs. |39],f|0|], we write down the system 
of equations for the time evolution of the pairwise correlation functions 



G AB (X,t) = - KC{t) [2G AB (\,t) + G AA {X,t) + G BB (X,t) 
+ 2 D A\ G AB (X, t) + I AB (X) + Tab (A), 



+ 



(13) 



9 



G AA (X,t) = - 2KC{t) G AB (X,t) + G AA (X,t) 



+ 



+ 2DA x G AA (X,t) + I5(X) + T AA (X), (14) 

G BB (X,t) = - 2KC(t) [G AB (X,t) + G BB (X,t)] + 

+ 2DA X G BB (X,t) + I 6(X) + T BB (A) , (15) 

where denotes the Laplace operator acting on the rf-dimensional variable A, the symbol 
Iab(X) in Eq.(13) describes the correlations in the production rates of A and B particles; it 
is equal to zero in the case I (uncorrelated generation of particles) and to 

Iab(X) = ld (X g )I8{\X\-X g ) (16) 

in case II, when particles are introduced into the system in correlated A — B pairs. Finally, 
in Eqs.(13) to (15) the terms T^- denote the correlation functions of the third order. 

The time evolution of the pairwise correlations is coupled to the evolution of the third- 
order correlations, which, in turn, depends on the correlations of the fourth order. Thus 
the non-linearity of the reaction-diffusion Eqs.(4) and (5) entails an infininite hierarchy of 
equations for correlation functions and one has to resort to some approximate methods. 

The most commonly used method invokes truncation of this hierarchy approximating the 
third-order correlation functions in terms of C(t) and Gij(X,t) f50| . Such an approach, as 



it was first noticed in Ref. |T4j], results in the Smoluchowski-type approximate results with 
improved numerical coefficients and is appropriate for the description of the single-species 
reactions A + A — >• O, but not for the description of reactions involving two different types 
of particles. The point is that such an approximation misses an important conservation 
law, which is specific for A + B — > O reactions. Namely, the reaction process conserves the 
difference Z(r,t) of local densities, Z(r,t) = C A (r,t) — C B (r,t), which changes only due to 
the diffusion of particles and thus is a pure diffusive mode of the system. Conservation of 
Z(r, t) entails, in turn, the conservation of the combination of pairwise correlation functions, 
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S-(X,t) = GAA(X,t) + (jbb(A, t) — 2(xab(A, £), which is also a pure diffusive modeQ. Con- 
sequently, only such truncation scheme should correctly describe the behavior of the binary 



reaction, which does not violate the important conservation laws [|14]]. 

The most simple truncation scheme, which preserves the conservation laws, has been 
first proposed in Ref. [IT|. In this scheme the third-order correlation functions, i.e. 7^, were 



set equal to zero. Such a truncation, as it was shown in Refs. 2!| and [38)-flQ is equivalent 



to the assumption that fields 5CA,B(X,t) have Gaussian distribution. Then, the fourth- 
order correlation functions automatically decouple into the product of pairwise correlation 
functions and the third-order correlations are equal to zero. Such an approach leads to, for 
instance, the correct long-time decay law of the particles mean densities, i.e. the t~ d ^ 4 -\aw, 
but fails to reproduce correctly the intermediate time behaviors; at intermediate times this 
approach predicts essentially the same behavior as the formal kinetic "law of mass action" 
and thus disregards the effects of particles' diffusion at the intermediate-time stage. 

This shortcomming has been revisited and improved in Refs. [P^ -iD|, in which it has been 



shown that correlation functions of the third order are small only in the limit |A| > R, while 
in the domain |A| ~ R they are singular and this singularity has an impact on the behavior 
of the pairwise correlation at the intermediate times. In a discrete-space picture, essential 
at scales |A| ~ R, the third-order correlation functions have been computed explicitly |39 . 



T AA « T BB « T AB ~ j d {R) C{t) S(\X\-R). (17) 

Substituting Eq.(17) into Eqs.(13) to (15) one gets then a closed with respect to C(t) 
and Gij(X,t) system of equations. To solve them, it is expedient to represent the pairwise 
correlations in the form 

Gij(X,t) = dij(\,t) + gij(X,t), (18) 

where Gij(X,t) denotes a "singular" part of the pairwise correlation functions, which ac- 
counts merely for the behavior of the third-order correlations, and gij(X,t) - the "fluc- 



2 One may readily verify [fj9) that Taa + Tbb — 2Tab is exactly equal to zero. 

11 



tuational" part, which accounts for the fluctuation spectrum of the external source and 
fluctuations stemming out of reaction and diffusive processes. 

The "singular" part of the pairwise correlation functions has been determined in Ref. 
f39|j . In particular, the leading at sufficiently large times behavior of Gjj(|A| = R, t) have 
been found to be as follows: 



Gij{,R,t) « C(t) (nt/SD) 1 / 2 , (19) 



in one-dimensional, 



<}^„C< t) !=M, (20) 



and 



Gij(R,t) « 6{t) 8irDR, (21) 

in two-, and three-dimensional systems, respectively. 

Now, inserting Eqs.(19) to (21) to Eq.(9) we obtain the following equation for the time 
evolution of particles mean density 

cw = - TT^k ^ w + to(|A| = *< () 1 + i + k'ikm (22 » 

where Ks(d) obeys, as t — > oo, 

K s (d = 1) « ( 8 -§y», (23) 

< 24 » 

and 

K s (d = 3) w STrDi? (25) 



One may readily notice that in three-dimensions the Kg{d), Eq.(25), coincides exactly 
with the so-called "diffusive" Smoluchowski constant; a reaction constant which has been 
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first calculated by Smoluchowski in his approximate description of the effects of diffusion 
on the chemical reactions kinetics. This constant accounts for, heuristically, the "resistivity" 
of random, diffusive transport of particles with respect to reaction. Employing the Smolu- 
chowski method, the analogues of such a constant have been obtained in Ref. |59| for one- 



and two-dimensional systems. Remarkably, our results in Eqs.(23) and (24) coincide with 
those obtained in Refs. U and [fj9| . We note also that the prefactor before the brackets in 



Eq.(22), i.e. the ratio K app = KKs{d)/(K + Ks(d)), is the so-called effective or apparent 
reaction constant, which was first derived for three-dimensional systems in Ref. [0]. There- 
fore, an account of the "singular" part of the third-order correlation function and subsequent 
extraction of the "singular" part in the pairwise correlators leads us to the results equivalent 
to those obtained in terms of the Smoluchowski approach. 

Hence, Eq.(12), in which one sets <7ab(|A| = R, t) = and K s (d) = oo reduces to the 
formal kinetic "law of mass action", while setting <7ab(|A| = R, t) = and using K$(d) 
as in Eqs.(23) to (25), one obtains the effective kinetic equation of the Smoluchowski- type 
approach. Below, we proceed to show that taking into account the time evolution of the 
pairwise correlations, i.e. the term ^abCA] = R, t), one arrives at completely different 
physical behavior as compared to the ones predicted by the formal kinetic and Smoluchowski 
approaches. 

Finally, we obtain the following system of equations, obeyed by the "fluctuational" part 
of the pairwise correlation functions. It reads 



9A B {\t) = - KC{t) 2g AB (\,t) + g AA {\,t) + g BB {X,t) 



+ 



+ 2DA x g AB {\t) + I AB (\), (26) 
g AA {\t) = - 2KC(t) \g AB {\t) + g AA {\t)] + 2 D A x g AA (\,t) + I 5(A), (27) 

9BB(\,t) = - 2KC{t) [g AB {X,t) + g BB {X,t)] + 2 D A A g BB {\, t) + J 5(A) (28) 

Equations (26) to (28), accompanied by Eq.(22), are now closed with respect to mean par- 
ticles' densities and pairwise correlations, and suffice the computation of the time evolution 

13 



of the monomer-monomer reaction scheme. 



III. KINETICS OF THE MONOMER-MONOMER REACTION SCHEME. 

Below we will analyse solutions of Eqs.(22) to (28) in systems of different dimensionalities 
and with different types of external particle generation. The derivation of results in case of 
one-dimensional systems will be presented in detail. The steps involved for such a derivation 
in higher dimensions are essentially the same and here we will merely discuss the results. 

A. Low dimensional systems. 

Let us start with the case of one-dimensional systems in which an external source pro- 
duces uncorrelated A and B particles. 

We note first that the system of equations (26) to (28) possesses two integrable combi- 
nations 

S_(\,t) = 2g AB (X,t) - g AA (X,t) - g BB (X,t), (29) 

which is related to the conserved property Z(r,t), and 

S+(X,t) = 2g AB (X,t) + g AA (X,t) + g BB (X,t) (30) 

These integrable combinations obey 

5_(A,t) = 2DA x S_(X,t) - 2I5{X), (31) 

which is thus the pure diffusive mode of the system, not affected by the reaction, and 

S+(X,t) = 2DA x S + (X,t) - AKC{t)S + {X,t) + 2 15(A) (32) 

The desired property, i.e. the correlation function g AB (X,t) which enters the Eq.(22), may 
be then expressed in terms of these integrable combinations as 

9A B {\t) = \(S.(X,t) + S + (X,t)) (33) 
14 



Consider now the solutions to Eqs.(31) and (32) in one-dimensional systems. Applying 
the Fourier transformation over the variable A, 

1 r°° 

S±(u,t) = —= \ dX exp(^A) S±(\,t) (34) 

to Eqs.(31) and (32), and assuming that at t = no A and B particles were present in the 
system, one readily gets that the Fourier-images of the integrable combinations follow 

S_(u,t) = - I J*dr exp(-2Druj 2 ) = 

= ~ n 2 J AT- ( X - exp(-2Dru; 2 )), (35) 
Du z y 2n 

and 



S + (uj,t) = - I J*dr exp(-2Druj 2 - AK J* dr' C{r')) (36) 

Now we notice that in the extreme situation, when reaction act occurs at any encounter of 
any A and B particle (i.e. when K = oo) the second integrable combination S + (u, t) vanishes 
since the integral f*dr' C(r') is obviously positively defined. One can show, however, that 
even for the finite Ks the influence of S + (u, t) on the pairwise correlation function, Eq.(33), 
is not essential at large times and the dominant contribution to g^B comes from S-(u,t). 

We note that setting t = oo in Eq.(35) we obtain that S^(u,oo) has a steady-state 
spectrum of the form S-(u, oo) ~ l/^ 2 , i.e. the spectrum which has an essential singularity 
when uj — > 0. Such a singular behavior of the fluctuation spectrum of the pairwise corre- 
lations in systems with binary reactions followed by an external uncorrelated production 
of the reactive species has been first obtained, in terms of a different than ours approach, 
in Refs. |L6j and |T?|]. The authors concluded thus that the steady state of such a system 



is highly anomalous; since such a singularity is not integrable in low dimensional systems, 
the steady state values of the integrable combination 5_(A, t) and thus of the correlation 
function g^B are infinitely large, which means that as time evolves the system progressively 
coarses into the domains containing particles of only one type. 

Consider now how the integrable combination S-(\,t) and the correlation function g^B 
grow in time. Taking the inverse Fourier transformation of the first line in Eq.(35) we get 
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s -^ i) ^-755i'$ exp( - AV8Br) (37) 

The integrand in Eq.(37) is a bell-shaped function with its maximum at point r = X 2 /8D. 
For bounded A, the bulk contribution to the integral comes from the algebraic tail 1/y/r 
and consequently, the leading at t ^> R 2 /8D behavior of the integrable combination follows 



S-(\\\ = R,t) = - (38) 

Accordingly, the absolute value of the "fluctuational" part of the pairwise correlation func- 
tion grows in time as 



^(|A| = J2,*) = -I^J^ (39) 

Inserting the just derived expression into the Eq.(22) and noticing that the correct asymp- 
totical behavior of the mean density obtains when the terms in brackets compensate each 
other, i.e., when 



C(t) = J-g AB (\X\ = R,t), (40) 



we find that 



C(t) = 1^ (41) 

i.e. in one-dimensional systems with random uncorrelated generation of the reactive species 
the mean particle density grows sublinearly in time as time progresses. The behavior as in 
Eq.(41) has been also obtained earlier in Refs. flTT|-f22]|. 



Consider now how the situation will be changed in the case II, when A and B particles 
are introduced into the reactive bath as correlated pairs. In this case one readily gets that 
the Fourier-image of the integrable combination S-(X,t) obeys the following equation 

S_(u,t) = -2Du 2 S-{u,t) - 2/(1 - cos(uX g )) (42) 

whose solution will read 
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= I (1 - cos^X,)) exp( _ 2DtQjl)) (43) 
Duo z y lit 

One may readily notice a very important feature of Eq.(43); in a striking contrast to the 

case I, here the steady-state spectrum is not singular in the limit to — > 0, but tends to a 
constant value 

J ^2 

S(u -> 0,t = oo) = =?-, (44) 



which means that S-(X, oo) and hence, #ab(A, oo) are bounded in systems of any dimen- 
sionality, and thus the well-defined steady-state mean density C(t = oo) also exists. We 
notice, however, that the steady-state pairwise correlation function is proportional to \ 2 g 
and thus may increase indefinitely with growth of X g . This unbounded growth is, of course, 
quite consistent with the result in Eq.(39), since the limit X g — > oo corresponds to the case 
of uncorrelated generation of particles. 

Now, the inverse Fourier transformation gives us 

m s I f l dr r . X 2 . 

s - A '* = " "7r^ / ~r l exp ~ 

\J2-kD Jo V r 8 -^ r 



1 (X-X g )\ 1 (A + A g ) 2 

= - (r(-i/2, 2L) _ Ir(-i/2,^^ 



ir ( -i/2,%±^)}, ,45) 



8L>t 

where T(a,x) denotes the incomplete Gamma function [[51 



Consider now the asymptotic behavior of the pairwise correlation function, Eq.(33), for 
different values of parameters A and X g , and different values of time t. 

One readily gets from Eq.(45) that at sufficiently short times, when A ^> X g 3> SDt, the 
pairwise correlation function obeys 

3IX 2 J8Dt) 3 / 2 A 2 , 

9ab(X, A 9 , t) = - 8V - mA «p(-— ), (46) 

which shows that correlations drop off as a Gaussian function at large scales. 
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Now, at short scales, such that A <C \ g and A <C 8Dt, and when A 9 is sufficiently large, 
\ g ^> 8Dt, we obtain that Pad (A, A 9 ,t) obeys Eq.(39), which is not a surprising result since 
at such scales the correlations in the particles' injection should be irrelevant. 

Within the opposite limit, when SDt <C A <C \ g the correlator follows 

^ x \ USDtf' 2 , A 2 . 

9AB {\Kt) ~ - exp(-— ), (47) 

which is reminiscent of the behavior in Eq.(46). 

Finally, in the limit when both A <C SDt and \ g <C 8Dt (and A 9 > A), i.e. in the limit 
of very long times, we find the following asymptotic expansion 

« - [1 - I -|=exp(- i A l) + 0(1/0], A s > A, (48) 

where the symbol 0(l/t) signifies that the correction terms decay with time as l/t. 

Equation (48) suffices to derive the large-t evolution of the mean particle density in the 
case of generation by correlated pairs, which reads 



^^^I'-i^w*^ 8 (49) 

Equations (49) reveals two surprising features; first, the steady-state density turns out 
to be dependent both on the diffusion constant and on the radius of pairs, generated by 
the source. Such a dependence is, of course, inconsistent with the predictions of the formal 
kinetic approach, based on the text-book "law of mass action". Second, the approach of 
particles' densities to their steady-state values obeys a power-law dependence, in a striking 
contrast to the exponential one predicted both by the formal kinetic and the Smoluchowski 
approach. 

To close this subsection let us briefly consider the behavior of solutions of the reaction- 
diffusion equations (22), (26) to (28) in two-dimensional systems. 

In the case I, we have from Eqs.(31) and (33) that as t — > oo the pairwise correlation 
function grows (by absolute value) as 

g AB {\\\ = R,t) w -Iln(Dt/R 2 ), Dt > R 2 , (50) 



and consequently, we get from Eq.(40) that in this case at large times the mean particle 



density exhibits logarithmically slow growth [16-22 



C(t) « \/lln(Dt/R 2 ) (51) 

Now, in the case II, we obtain that the steady-state exists and <7ab(|A| = R, oo) behaves as 

g AB (\\\= R,oo) « -Iln{\ g ), (52) 

and hence, the steady-state density turns to be a slowly growing function of the radius of 
the generated pairs, 



C(t = oo) « yjl ln(X g ) (53) 
Some analysis shows also that such a steady-state is approached via an algebraic law 

C(t) - C(t = oo) « (£>i)~\ (54) 
in contrast to the exponential in time dependence predicted by mean-field descriptions. 

B. Three-dimensional systems. 

As we have already mentioned, in the case I the steady-state fluctuation spectrum is 
characterized by an essential singularity of the type 1/u 2 as uj —>■ 0. In three-dimensional 
systems such a singularity is integrable, which insures that the steady-state correlations 
exist and vanish as A — > oo. Consequently, the steady-state mean particle density C(t — oo) 
should exist in 3D also in this case. Let us analyse now the form of these correlations. 
Solving Eq.(31) in 3D we find (up to the correction terms which are exponentially small 
with K) that 

g AB {\t = oo) w - g^X> ( 55 ) 

i.e. A — B correlations vanish with the distance between particles A as 1/A, which shows that 
in the monomer-monomer catalytic scheme taking place in three dimensional systems the 
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correlations in the steady-state show a quasi-long-range order decaying as the first inverse 
power of the interparticle distance. 

Now, substituting Eq.(55) into the Eq.(22) we find the following expression for the steady- 
state density in 3D, 



= = V<£Z« + (56) 
which shows that A — B correlations lead here to an effective renormalization of the reaction 



constant in the steady-state, i.e. C(t = oo) has the form C(t = oo) = J I /K app , where K app 
is the mentioned above apparent reaction constant 0. 

Consider now how such a steady-state is approached at long times. Expanding the 
solution of Eq.(31) near the steady-state, we have that pairwise correlations approach the 
steady-state as a power-law, 

g AB (\M=R,t) « — ^ [l " -4= + 0(l/t)}, (57) 
which yields, in turn, a power-law relaxation of the mean particle density to the steady-state 

C(t) - C(t = oo) « (Dty 1/2 (58) 

Therefore, in contrast to low dimensional systems, in three dimensional systems with 
random uncorrelated generation of the reactive species the essential singularity in the fluc- 
tuation spectrum is integrable, correlations vanish with the distance between particles and 
the steady-state mean particle density exists. However, the steady-state density is differ- 
ent from that predicted by the mean- field "law of mass action" and shows, in particular, 
dependence on the particles' diffusivity D. Besides, Eq.(58) reveals that here approach to 
the steady-state is described by a power-law with the characteristic exponent —1/2, which 
is essentially non-mean-field behavior. 

Finally, for the case II we find the following results for the correlation function and mean 
density. In the steady-state the A — B correlations are equal to zero for A > X g (again, apart 
from the exponentially small with K terms). In the domain A < X g the correlations exist 
and are described by 
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g AB (X,t = oo) » - (1 - (59) 

which reduces to the result in Eq.(55) when \ g = oo. In constrast to the behavior as in 
Eq.(55), however, the correlations vanish at finite values of the correlation parameter A. 
Now, Eq.(59) yields for the steady-state mean-particle density 



c < = »> = f 4 + S^'- < 60) 

which is less than the steady-state density in the case I, Eq.(56), due to a factor 1 — R/\ g , 
which renormalizes the Smoluchowski constant. Consequently, for \ g > R apparent rate 
constant here takes the form 

_ 8ttDRK 
app ~ 8irDR+(l-R/\ g )K (6 } 

We find also that such a steady-state is approached via a power-law, 

C(t) - C(t = oo) « (Dt)-^ 2 (62) 

which is faster than the approach described in Eq.(58), but still very different from the 
exponential behavior predicted by mean-field analysis. 



IV. CONCLUSION. 

To summarize, we have shown that both the cases I and II the fluctuations effects dom- 
inate the kinetics of the monomer-monomer catalytic scheme involving diffusive particles 
and induce essential departures from the predictions of the mean-field approaches. In the 
case I, the effects of fluctuation are especially pronounced in low dimensional systems - the 
steady-state does not exist and mean particle density grow indefinitely in time, in absence 
of hard-core exclusion between particles. In three dimensions the steady-state exists, but is 
characterized by very strong interparticle correlations, which, in turn, have a strong impact 
on the value of the steady-state mean particle density. The steady-state density is different 
from that predicted by mean- field "law of mass action". The approach to this steady-state 
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is described by an anomalous power-law with the characteristic exponent —1/2, which stems 
from the presence of essential singularity in the steady-state fluctuation spectrum. In the 
case II, the steady-state fluctuation spectrum and the steady-state mean particle density ex- 
ist in any dimension, but show an anomalous, non-mean-fields dependence on the particles' 
diffusivity and the radius of pairs, generated by the source. Approach to the steady-state 
follows a universal power-law with the characteristic exponent —d/2, which resembles, apart 
from the dependence of the prefactors on the system parameters (e.g. constant of the back- 
ward reaction) the long-time approach to the equilibrium in reversible chemical reactions 



12ILM38H421. The origin of such a behavior is that the fluctuation spectrum in the steady- 



state is flat at small values of the wave- vector, i.e. the essential singularity in the steady-state 
spectrum of fluctuations is screened. 
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Figure Caption. 

Fig.l. Langmuir-Hinshelwood reaction on a two-dimensional catalytic substrate. Black 
and grey spheres denote particles of A and B species, respectively; (1) describes the situation 
in which an A and a B appear within the reactive distance from each other and may enter 
into reaction. 
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Fig. 1. Argyrakis et al., J. Chem. Phys. 



